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ABSTRACT 

A semi-analytic model is proposed that couples the Press-Schechter formalism for the 
number of galaxies with a prescription for galaxy-galaxy interactions that enables to 
follow the evolution of galaxy morphologies along the Hubble sequence. Within this 
framework, we calculate the chemo-spectrophotometric evolution of galaxies to obtain 
spectral energy distributions. We find that such an approach is very successful in 
reproducing the statistical properties of galaxies as well as their time evolution. We are 
able to make predictions as a function of galaxy type: for clarity, we restrict ourselves 
to two categories of galaxies: early and late types that are identified with ellipticals 
and disks. In our model, irregulars are simply an early stage of galaxy formation. In 
particular, we obtain good matches for the galaxy counts and redshift distributions 
of sources from UV to submm wavelengths. We also reproduce the observed cosmic 
star formation history and the diffuse background radiation, and make predictions as 
to the epoch and wavelength at which the dust-shrouded star formation of spheroids 
begins to dominate over the star formation that occurs more quiescently in disks. A 
new prediction of our model is a rise in the FIR luminosity density with increasing 
redshift, peaking at about z ~ 3, and with a ratio to the local luminosity density 
Pl,u{z — Zp ea k)/ Pl,u{z = 0) about 10 times higher than that in the blue (B-band) 
which peaks near z ~ 2. 
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1 INTRODUCTION 

In contrast with the local universe, where only 30 % of the 
bolometric luminosity is released in the IR/submm wave- 
length range (Soifer & Neugebauer 1991), there is a grow- 
ing amount of evidence that the high-redshift universe was 
much more opaque. Indeed, the discovery of the Cosmic In- 
frared Background (CIRB) at a level ten times higher than 
the no-evolution predictions based on the IRAS local IR lu- 
minosity function, and twice as high as the Cosmic Optical 
Background obtained from optical counts, has shown that 
dust extinction and emission are key processes for high- 
redshift galaxies (Puget et al. 1996; Guiderdoni et al. 1997; 
Fixsen et al. 1998; Hauser et al. 1998; Schlegel, Finkbeiner 
& Davis 1998). Deep surveys with the ISO satellite at 15 (im 
(Oliver et al. 1997; Aussel et al. 1999; Elbaz et al. 1999) and 
175 fim (Kawara et al. 1998; Puget et al. 1999), and with the 
SCUBA instrument at 850 /im (Smail, Ivison & Blain 1997; 
Barger et al. 1998; Hughes et al. 1998; Eales et al. 1999; 
Barger, Cowie & Saunders 1999) have begun to resolve the 
CIRB into its brightest contributors. Although identification 
and spectroscopic follow-up of submm sources are not easy, 



such studies seem to reach the conclusion that an important 
fraction of these sources are the high-redshift counterparts 
of the local luminous and ultraluminous IR galaxies (LIRGs 
and ULIRGs) discovered by IRAS (Smail et al. 1998; Lilly et 
al. 1999; Barger et al. 1999). In the optical and near-IR win- 
dow, careful examinations of the Canada-France Redshift 
Survey (CFRS) galaxies at z ~ 1, and Lyman break galaxies 
at z ~ 3 and 4 have revealed a significant amount of extinc- 
tion (Flores et al. 1999; Steidel et al. 1999; Meurer, Heckman 
& Calzetti 1999). This has lead to a reassessment of previ- 
ous estimates of the UV fluxes, and consequently of the star 
formation rates in these objects, which are now found to 
be higher by a factor 2 to 5. In light of these observations, 
one can view the far infrared background as a sink for the 
hidden aspects of galaxy formation. At optical wavelengths, 
ellipticals and spheroids are old, even at z ~ 1. No evidence 
is seen for either the luminous formation phase, or the early 
evolution at these wavelengths. One therefore draws the con- 
clusion that ellipticals and, more generally, most spheroids 
must have formed in dust-shrouded starbursts. 

In this paper, we show how separate tracking of disk and 
spheroid star formation within a simple but explicit cosmo- 
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Table 1. Parameters of the cosmological model, h is the Hubble 
constant at present in units of 100 km/s/Mpc. ag is the rms mass 
fluctuation in spheres of Sh~^ Mpc. 



n 


Ao 


h 


0"8 


n B 


0.3 


0.7 


0.65 


0.86 


0.02 h~ 2 



logical framework enables one to substantiate this view and 
show that current data is consistent with spheroid forma- 
tion dominating the submillimetre background. In section 
2, we briefly describe how we connect the various pieces of 
our semi-analytic model together, to differentiate between 
spheroids and disk galaxies and use the STARDUST model 
galaxy spectra (Devriendt, Guiderdoni & Sadat 1999) to 
make the link with observations. The approach is very simi- 
lar in spirit to that of Devriendt & Guiderdoni (2000) - DG00 
hereafter, so we point out the differences/improvements with 
respect to this previous work. Section 3 presents the results 
obtained with this new implementation of collision-induced 
galaxy formation. Finally, we put our work in context with 
other studies and draw conclusions in Section 4. 



2 OVERVIEW OF THE MODEL 
2.1 Cosmological framework 

To match the combined observations of the cosmic mi- 
crowave background anisotropies on sub-degree scales by 
BOOMERanG (de Bernardis et al. 2000) and MAXIMA 
(Balbi et al. 2000), and of high-redshift supernovae (Riess 
ct al. 1998; Schmidt et al. 1998; Perlmutter et al. 1999), 
we have chosen a flat universe model with a cosmological 
constant. More specifically, the ratio of total matter density 
to critical density and the reduced cosmological constant 
are taken to be fio = 0.3 and Ao = 0.7 respectively. The 
Hubble constant Ho was chosen to lie within the error bar 
quoted by the Hubble Space Telescope key project value 
(Freedman et al. 2001). The value of n B = 0.02 h~ 2 is the 
one currently favored by various cosmological probes of Big 
Bang nucleosynthesis (Bania et al. 2002). Finally the nor- 
malization of our non-tilted primordial Harrison-Zeldovich 
power spectrum as = 0.86 ensures fair agreement both with 
the amplitude of the cosmic background radiation measured 
by COBE and the local cluster abundance. We note that 
recent cluster surveys seem to favor slightly lower values 
around ug ~ 0.73 (e.g. Lahav et al. 2002) but our conclu- 
sions are quite insensitive to a change of 20 % in any of the 
parameters presented in table 1. 

In the semi-analytical ab initio approach that we out- 
line in the next section, galaxies form from Gaussian ran- 
dom density fluctuations in the primordial matter distribu- 
tion, dominated by cold dark matter (CDM). To follow a 
bound perturbation, we use the well known spherical top- 
hat collapse model: a homogeneous spherical perturbation 
grows along with the expanding universe, until self-gravity 
results in turn around and (non-dissipative) collapse. We 
then compute the number of perturbations of a given mass 
which collapse at a given redshift using the standard Press 
& Schechter (1974) prescription and starting with a primor- 
dial power spectrum for the density field given by Bardeen et 



al. (1986), corrected for the presence of baryons (Sugiyama 
1995). To take into account mass build up of halos through 
mergers, we employ the extended Press-Schechter formalism 
described in e.g. Lacey & Cole (1993), in a similar fashion as 
in Balland, Silk & Schaeffer (1998) - BSS98 hereafter. Note 
that this is slightly different from Devriendt & Guiderdoni 
(2000), where the prescription to compute the number of col- 
lapsed haloes was given by the peaks formalism developed 
in Bardeen et al. (1986) and where halo merging histories 
were not considered. 



2.2 A simple semi-analytic model 

To make contact with observations, star formation and spec- 
tral evolution must be incorporated. We follow the general 
framework described by White & Frenk (1991), that is to 
say, within virialized dark matter halos, we let gas cool ra- 
diatively, settle into a disk and form stars. Technical de- 
tails about the semi-empirical recipes used to model the dif- 
ferent astrophysical processes are given in DG00. Chemical 
and spectral evolution are then computed with the STAR- 
DUST stellar population synthesis model from Devriendt et 
al. (1999). Detailed references can be looked up in this paper, 
but, as our results in the far-IR/submm crucially depend 
upon it, we summarize the basic assumptions/ingredients 
concerning dust that go into this modelling. The first of these 
is that metals and dust are homogeneously mixed with stars 
in galaxies of oblate ellipsoidal shape. The extinction curve 
of each galaxy is then assumed to scale with its gas and metal 
contents, in a manner that is calibrated on the Milky Way 
and nearby spirals. From such a model, the radiative transfer 
can be computed analytically, and yields the total amount 
of stellar light absorbed at UV/optical wavelength. Follow- 
ing energy conservation, this energy must be re-emitted 
at IR wavelengths. Implementing the 3-component absorp- 
tion/emission model of Desert, Boulanger & Puget (1990), 
(a mixture of PAHs, very small grains and large grains) as 
in STARDUST, i.e. in combination with the empirical correla- 
tion of the bolometric IR luminosity with IRAS colours, we 
finally link the optical and the far-IR/submm windows in a 
self-consistent way. 

We emphasize that, as discussed in detail by DG00, 
there are three key parameters in such models, the star for- 
mation efficiency the feedback efficiency €sn and the 
extent of the gaseous disks / c . We take parameter values 
that are fairly similar to the ones used by these authors, 
with ft = 60, esN ~ 0.2 and f c — 5. These values are well 
within the uncertainties of observations by Kennicutt (1998) 
and Bosma (1981) for the star formation efficiency and the 
gaseous extent of cold disks, respectively. Numerical simula- 
tions by e.g. Thornton et al. (1998) tend to give values closer 
to 0.1 for feedback efficiency; however our higher value is a 
consequence of an attempt to reduce the number of small 
objects overestimated by the Press-Schechter prescription, 
as well as a means to cure the "cooling catastrophe" , where 
too much gas cools in low mass halos at high redshift. The 
qualitative effect of each of these parameters can be summed 
up in the following way: 

• Increasing f3 decreases the normalisation and faint-end 
slope of the optical and IR counts, because star formation 
is reduced and takes place at lower redshift. 
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• Increasing esiv decreases the normalisation and faint- 
end slope of the optical and IR counts by quenching star 
formation in higher and higher mass galaxies that form on 
average at lower and lower redshifts. 

• Increasing f c means increasing the normalisation of the 
optical counts, and decreasing that of the IR counts, since 
this is equivalent to reducing extinction. 

The main drawback of DGOO is their modelling of ul- 
traluminous dusty starbursts: they use an ad hoc fraction of 
massive objects F, calibrated locally on the IRAS luminosity 
function and which is evolving with redshift as F oc (1 + z) 6 . 
In this paper, we replace this "recipe" by a physically mo- 
tivated collision model which allows us to differentiate be- 
tween disks and ellipticals. 



2.3 Modelling galaxy types 

To separate disk galaxies from spheroids we use the prescrip- 
tion proposed in BSS98. It is based on the amount of energy 
exchanged by a galaxy with its neighbours via gravitational 
interactions during its lifetime. Interactions are modelled as 
rapid non-merging collisions following the Spitzer impulse 
approximation (Spitzer 1958). If we denote by A the in- 
stantaneous relative rate of energy exchange, (l/E)dE/dt, 
at epoch t, the cumulative effect of interactions occurring 
throughout the galaxy lifetime (from its birth at the red- 
shift Zni at which its host halo collapses to the redshift z of 
observation) is obtained by integration: 



k dt , 
A—dz 

dz 



(1) 



In the above expression, the main parameters describing 
an interaction (impact parameter, relative velocity between 
the interacting galaxies, relative mass, number density of 
neighbouring galaxies) enter through the quantity A. The 
scalings with these parameters are derived from the numer- 
ical simulations of Aguilar & White (1985). In this way, it is 
possible to overcome the difficulty of modelling analytically 
the details of the interaction, which is a highly non-linear 
process. 

The effects of cosmology enter equation (1) via the con- 
version of the loockback time t into redshift z and the scaling 
of the interaction parameters with redshift. For instance, in 
the present cosmology (see table 1), the relative velocity be- 
tween two interacting galaxies is given by: 



vo 



^n (i + z) + \ (i + z)- 2 . 



(2) 



The evaluation of equation (1) for a galaxy in the field 
(see BSS98 for details) gives (in a flat, fio+Ao = 1, universe): 



A{z) = A*n { 



(l + z»i) 3 -(1 + z) 3 

n 2 
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fi 2 
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Znl 



(fio(l + z n i)' A + A )(n (l + z) 3 + A ) 



with A* 
that in the limit Ao 



»0(1 + Zni) 3 + A l "I 

n (l + 2) 3 + A„ J J ( ' 

a 1.35 x 1CT 4 for h = 0.65. It is easy to verify 



A*[(l + z ni f - (1 + z) 3 ] which is the result for an fi () = 1 
universe (BSS98). 

Morphological types are then defined according to the 
value of A using the following rules: 

• Galaxies for which A is lower than a certain threshold 
value A sp i are identified as disks: they correspond to galaxies 
having experienced few, if any, interactions and consequently 
have a small A. 

• Galaxies for which A is higher than a certain threshold 
value A e ii > A sp i are identified as ellipticals: they corre- 
spond to galaxies having undergone substantial energy ex- 
change during their lifetime. 

• Galaxies such that A spi < A(z) < A e « are identified 
as SO galaxies. 

The values of A sp i and A e u are fixed by requiring that 
the model produces the observed amount of each type in the 
field at the present epoch. Typical values are A sp i « 3.10 -3 
and A e ii ~ 10~ , and are not very sensitive to the cos- 
mological model assumed (standard or ACDM). The above 
conditions on A can then be inverted to obtain conditions on 
the formation redshift z„i(z) of galaxy types. It is important 
to realize that the formation redshift of a galaxy of a given 
type at epoch z depends on z. That is to say that in our 
model the morphology of a galaxy is not fixed once and for 
all by initial conditions but, on the contrary, it evolves pro- 
gressively from late-type to early-type. However, contrary 
to most other semi-analytic approaches, this progression is 
monotonic: once a spiral galaxy has become SO it can only 
evolve to the elliptical stage, and it can only remain ellipti- 
cal once it reaches that stage. The condition for a galaxy to 
be elliptical at redshift z is that its formation redshift z n i is 
such that: 



Z„l > Z e ll(z) 



(4) 



where z e u(z) is defined by equating A(z) to A e u and solving 
for z n i. Equivalently, spirals at epoch z are those galaxies 
that turned non-linear at an epoch such that: 



Znl ^ Zspii^Z) 



(5) 



where z spi (z) is defined by equating A(z) to A api . Finally, 
SO galaxies are such that: 



,i(z) < Znl < Z e ll(z). 



(6) 



equation (3) reduces to A(z) = 



Figure 1 illustrates the way morphological types are de- 
fined in this model. The formation redshift is shown as a 
function of the epoch z considered. The solid curves are ob- 
tained for A(z) = A e ii (upper) and A(z) = A spi (lower). El- 
lipticals occupy the upper region of the (z n i, z) plane, while 
disks occupy the lower region. SOs are located in between. 
The hatched area corresponds to z n i < z and is irrelevant. 

Consider the example of a galaxy that became non- 
linear at, say, z n i = 3. Figure 1 shows that up to redshift 
z ~ 2.8 this galaxy is identified as a disk. Between z ~ 2.2 
and z ~ 2.8, the same galaxy has experienced sufficient en- 
ergy exchange with its neighbours to be identified as an SO 
galaxy, but not quite enough to be an elliptical. By z ~ 2.2, 
it has become an elliptical. 

Depending on the previously identified morphological 
type, we decide which galaxies will undergo an "obscure 
starburst" . In other words, this phase is triggered when the 
energy imparted by encounters/collisions is strong enough to 
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12 3 



z 

Figure 1. (z n i,z) plane. The upper solid curve is for A = A e ;; 
and the lower one for A = A sp ;. The upper region is the locus for 
galaxies having experienced high energy exchanges due to gravi- 
tational interactions since they turned non-linear (the ellipticals 
region according to the present model), the lower region is the 
locus for low energy exchange (the disks region). In between, one 
finds the region for SOs. The hatched area is the region for which 
z > z n i and is therefore not populated by galaxies. 

alter the morphology for the first time, i.e. when late-type 
(spiral) galaxies first become early type (SOs). The transition 
from SO to elliptical is then supposed to happen smoothly, 
without another burst of star formation, as the gas poor SO 
galaxies are assumed to have reheated their surroundings 
enough during the burst to not be able to replenish their 
gas supplies through accretion of new material. The inten- 
sity and duration of this LIRG/ULIRG phases, are therefore 
controlled by the amount of gas available for star formation 
and the size of the galaxy at the time when the morphol- 
ogy change occurs. In practice, this dark phase is modelled 
by setting our three key parameters to /3 = 1, / c = 1 and 
€s n = 0.5, corresponding to high star formation efficiency, 
high dust opacity, and high feedback efficiency respectively. 
It is therefore completely coupled to the starburst, and as 
the luminosity of this latter decreases for lack of fuel (gas) 
supply, the optical depth of the galaxy goes down as well, 
revealing more and more of the stellar population at optical 
wavelengths. 



(~ 90 %) of the SCUBA sources are classified as early types 
in our model. The reason for such a change of behaviour lies 
in the well known negative k-correction, which makes galax- 
ies of the same bolometric luminosity as bright for the ob- 
server at redshift 5 as at redshift 0.5. This is only important 
in the submm (here for SCUBA at 850 microns), because 
the peak emissivity of dust in the source rest frame is be- 
tween 60 and 100 microns. Therefore, as our SOs/ellipticals 
approximately form at z > 2, the corresponding maxima of 
emission must be redshifted to wavelengths greater than 180 
and 300 microns respectively. 

This result is quite robust in the sense that its quali- 
tative features do not depend on the cosmological parame- 
ters. However, quantitatively, there is a marked difference: 
the domination of late type galaxies is more marked in a 
SCDM model, where galaxies tend to form later on average 
(see Silk & Devriendt 2000). This remark also applies to 
the model of DG00, where the phenomenologically evolving 
ULIRG fraction was the dominant contributor in the far-IR 
(175 microns). 

Although the general agreement of our predicted counts 
with the multi-wavelength data seems quite impressive, 
there are several caveats. At 15 microns, for example, one 
would say that we match the integral counts fairly well (up- 
per left panel of figure 3). But looking more closely, we can- 
not reproduce the change of slope seen in the ISOCAM dif- 
ferential counts (fig. 4). There are at least a couple of rea- 
sons why this could happen. First, the SEDs of the ISOCAM 
galaxies are different in the mid-IR from the ones used here 
as a template, which are based on IRAS observations of the 
local universe. Differential counts are very sensitive to the 
exact shape of the PAH features, which are crudely mod- 
elled here. One could also imagine that the discrepancy is 
partly due to the grain size distribution/chemical composi- 
tion evolution with redshift, as the redshift distribution of 
ISOCAM galaxies has a median z ~ 0.7. We plan to in- 
vestigate these issues in more detail but that is beyond the 
scope of this paper. Finally, because of the way interactions 
are modelled, each early type galaxy undergoes a starburst 
after its host halo has just collapsed, and it is not obvi- 
ous that ISOCAM sources (see figure 3) in which the vast 
majority are LIRGs (not ULIRGs), are properly described 
by such a violent process. Dynamical interactions (which 
are not modeled in detail here), triggering multiple milder 
starbursts, with time delays between them, might provide a 
more realistic description of these sources and this is another 
issue we plan to investigate in the future. 



3 RESULTS 

3.1 Galaxy counts 

Looking at figures 2 and 3, one realises that (except in 
the near IR bands), late type galaxies dominate over early 
types (compare the dotted curves with the dashed ones in 
each panel). This domination extends down to the far-IR, 
with the late-type galaxies still dominating the 175 micron 
ISOPHOT counts. At longer wavelengths, however, there 
is a dramatic change: the early- type galaxies completely 
swamp the contribution from late types. Indeed, one can see 
on the bottom right panel of figure 3 that the vast majority 



3.2 Redshift distributions 

Another constraint on our models comes from the redshift 
distributions. Their shapes seem to quite nicely match the 
observations in the I band with a mean redshift of the dis- 
tribution of ~ 0.6, (see top left panel of figure 5). In the 
far-IR (60 microns), the agreement with data gathered in 
the north ecliptic pole region (NEPR) is also fairly convinc- 
ing. As predicted in Silk & Devriendt (2000), inclusion of 
the cosmological constant, A, has shifted the near-IR and 
60 micron peaks towards higher redshifts and produced a 
high-redshift tail in the I band, bringing the models into 
closer agreement with the data. Disks and early-types are 
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Figure 2. UV/Optical/Near-IR counts from our analytic model, compared to existing data. The dotted curves stand for late type 
galaxies (disks); the dashed ones for early type galaxies (spheroids); and the solid ones for the sum of both contributions. Data is from 
Hogg et al. (1997) (U band), Williams et al. (1996) (F300W A B, B & I bands), Arnouts et al. (1997) (B band), Bertin & Dennefeld (1997) 
(B band), Gardner et al. (1996) (B, I & K bands), Metcalfe et al. (1996) (B band), Weir et al. (1995) (B band), Smail et al. (1995) (I 
band), Crampton et al. (1995) (I band), Moustakas et al. (1997) (K band), and Djogorvski et al. (1995) (K band). 



found in comparable proportion in the I-band, with a slight 
domination of disks for z > 0.3 and up to z = 1.5. 

There are at least a couple of reasons for this behaviour. 
First, in this redshift range, spheroids are already old, their 
star formation rates are very low and therefore their I-band 
luminosity comes from an old and dim stellar population. 
Secondly, spheroids in a massive starburst phase at these 
redshift experience high dust absorption, which reduces sig- 
nificantly their luminosity even in the I band. Going deeper 
in I magnitude leads to the selection of higher redshift galax- 
ies and we expect that the starburst population to eventually 
dominate the distribution around z « 2. 

We draw the reader's attention to the dramatic change 
occuring in the redshift distribution of sources when go- 
ing from ISOPHOT to SCUBA, i.e. from the far-IR, to the 
submm window (bottom panels of figure 5). In the far-IR, 
the vast majority of sources is predicted to lie at fairly low 
redshifts, whereas in the submm, there is a characteristic 
double peaked distribution. The first small peak at low red- 
shift is mainly due to disks, as can be seen on the lower 
right panel of figure 3 (A e // = 850 /im). Indeed, the Eu- 
clidian tail of the total counts (figure 3) at fluxes above 



~ 10" 1 Jy on this figure is clearly due to the domination 
of local bright disks (dotted line). In contrast, one notes on 
the lower-right panel of figure 5 an overwhelming domina- 
tion of the (fainter) sources located at high redshifts, as is 
also clear from figure 3. Once again, the main driver for this 
effect is the negative k-correction boosting observed fluxes in 
the submm wave-bands. Contrary to Devriendt and Guider- 
doni (2000) and Lagache et al. (2002), the second peak in the 
redshift distribution of 175 microns sources (around z = 1, 
see e.g. Lagache et al. (2002)) is missing from the model dis- 
cussed here. This could either be a result of the STARDUST 
SEDs being too warm for LIRGs, or of a too rapid evolu- 
tion in the number of LIRGs due to our simple modelling of 
galaxy interactions. Here too, triggering several starbursts 
per galaxy with delays between them, as happens in realistic 
galaxy interactions, could alleviate the problem. 



3.3 Diffuse background radiation 

The background light, as integrated from the multi- 
wavelength counts, seems to closely match the detection 
by Puget et al. (1996) and Hauser et al. (1998) in the 
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Figure 3. Infrared/submm counts from the analytic model described in the text. Coding for the lines is dotted: late type galaxies 
(spirals); dashed: early type (ellipticals, SOs); and solid: sum of both contributions. Data are from Elbaz et al. (1999) (15 fira), Kawara 
et al. (1998), Puget et al. (1999) and Dole et al.(2001) (175 (an), Smail et al. (1997), Eales et al. (1999), Barger et al. (1999), Scott et 
al. (2002) (850/im). 



far-IR/submm window (see figure 6). This tells us that we 
are not grossly overestimating or underestimating the faint 
counts, and that the global galaxy luminosity budget from 
the UV to the submm is likely to be correctly computed. An 
interesting remark is that this plot clearly shows that one 
has to go to the submm, around 1500 GHz (about 200 mi- 
crons), to see the early- type galaxies dominate over the late 
types. This is a factor 2 shorter in wavelength than found for 
the SCDM model. At any wavelength shorter than this (ex- 
cept in the near-IR between 2-5 microns), quiescent galaxies, 
forming up to ~ 10 Mq of stars per year dominate the total 
light emission (as well as the counts). Again that was not 
the case with the SCDM model where spirals dominated the 
whole background below 400 microns. Finally, we point out 
that the model marginally underestimates the background 
measured by DIRBE. This discrepancy could stem from a 
fraction of the population of galaxies having slightly differ- 
ent spectral shapes in the far-IR than our model assumes 
they have, or having a redshift distribution different from 
the one we predict they have. Of course, as these effects are 
degenerate, in the sense that colder sources will emit more 
flux in the far-IR than hot ones, as will closer sources, it is 



very difficult to decide which is the dominant without fur- 
ther guidance from the observations. 



3.4 Multi-wavelength comoving luminosity 
densities 

Figure 7 shows the comoving luminosity densities predicted 
by the model from the far-UV (bottom curves) to the far-IR 
(top curves). It is clear from the figure that there already 
is a greater amount of evolution in the far-IR than in the 
UV-optical, and we know from the counts (fig. 3) that the 
submm must show still more evolution than the far-IR. The 
agreement with the data is quite satisfactory both at low 
and high redshift and for all wavelengths. This tells us that 
the luminosity budget of our galaxies is relatively accurate. 
Another striking feature is that the model predicts that the 
spheroid (whether in dusty starburst phase or not) contri- 
bution to the comoving luminosity density becomes approx- 
imately equal to that of late type galaxies around redshift 2 
in the UV and far-IR and as early ar z ~ 1 in the B-band. As 
a result, the peak of the comoving luminosity density shows 
a complex behavior as a function of wavelength, starting 
around z ~ 3 in the UV and falling to z ~ 2 in the optical 



© 0000 RAS, MNRAS 000, 000-000 



Collision-induced galaxy formation: semi- analytical model and multi-wavelength predictions 7 




Figure 4. Infrarcd/submm differential counts from the model. Coding for the lines is dotted: late type galaxies (spirals); dashed: early 
type (ellipticals, SOs); and solid: sum of both contributions. Data are taken from various ISO surveys described in Elbaz et al. (1999) at 
15 fim, and various analysis of IRAS samples (see e.g. Guiderdoni et al. 1998 at 60 ^m). 




Figure 5. Panchromatic redshift distributions of galaxies predicted by our analytic model. Dotted curves represent late type galaxies 
(spirals); dashed: early type galaxies (ellipticals, SOs); and solid: the sum of both contributions. Data in the I band (top left panel) is 
from Crampton et al. (1995) (Canada-France Redshift Survey), whereas at 60 microns, it comes from Ashby et al. (1996) (North Ecliptic 
Pole Region). In both cases, the model curves have been renormalized to the total number of observed galaxies. 
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Figure 6. Diffuse background light emitted by star forming 
galaxies, from the UV to the submm. Coding for the lines is 
dotted: late type galaxies (disks); dashed: early type galaxies 
(spheroids); and solid: sum of both contributions. The data 
sources are indicated on the figure. 
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Figure 7. Comoving luminosity densities as a function of red- 
shift. As usual, coding for the lines is dotted: late type galaxies 
(disks); dashed: early type galaxies (spheroids); and solid: sum of 
both contributions. The three sets of curves represent predictions 
of our model at 2800 A(bottom triplet), 4400 A(middle triplet) 
and 60 microns (top triplet). Data sources are indicated on the 
figure. 
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Figure 8. Comoving star formation rate as a function of red- 
shift. The different curves represent the contribution of late type 
galaxies (spirals: dotted line); early type galaxies (ellipticals, SOs: 
dashed line); and the sum of both contributions (solid line). 
Empty symbols give raw measurements of different teams (renor- 
malizcd to match the cosmological model defined in table 1): cir- 
cles (Madau, Pozetti & Dickinson 1998); squares (Steidel et al. 
1999); triangles (Lilly et al. 1996) and pentagons (Connolly et al. 
1997). Filled diamonds are data from Gronwall (1998); Tresse & 
Maddox (1998); Treyer et al. (1998); Cowie, Songaila & Barger 
(1999); Steidel et al. (1999) and are dust-corrected as described 
in Somerville et al. (2001). The point at z ~ 0.7 is derived from 
far-IR ISOCAM sources (Flores et al. 1999). 

near-IR, before shifting again to higher redshifts (z ~ 3) in 
the far-IR and submm. This reflects the fact that both UV 
and far-IR luminosities are extremely sensitive to the instan- 
taneous SFR, whereas the optical and near-IR luminosities 
are tainted by a non-negligible contribution coming from an 
old stellar population. 

3.5 Cosmic star formation history 

The predicted comoving star formation rate is compared 
to observations on figure 8. Empty symbols are observed 
estimates at various redshifts, uncorrected for dust extinc- 
tion. Filled symbols indicate dust corrected estimates follow- 
ing the prescription of Somerville, Primack & Faber (2001) 
based on a luminosity dependant correction likely to be more 
realistic than the traditional scaling by a constant factor (be- 
tween 3 to 5) at all z (see their figure 9). The filled point 
at z ~ 0.7 is derived from infrared ISOCAM observations 
(Flores et al. 1999). 

The comoving SFR is in fair agreement with the mea- 
sured one, and the reasonable amount of K-band light pro- 
duced (figure 2) makes us confident that we do not over- 
produce either stars or metals. However note that at low 
redshift, the agreement is only for non-dust-corrected data, 
so the contribution of our spirals (which closely match the 
non corrected Madau curve at all redshifts) might be un- 
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derestimated in our model. From the perspective of metals 
only, the star formation rate could be higher, provided that 
a fair fraction of the heavy elements is ejected into the IGM. 
The model however predicts a peak in star formation around 
redshift ~ 3.5 (this is also true for the SCDM case), sensi- 
bly higher than what is usually assumed, current data being 
compatible with a flat star formation rate at redshifts z > 
1. A more detailed comparison with the SCDM cosmology 
used in Silk & Devriendt (2000) reveals that the comoving 
star formation rate increase from z — to z w 3 is steeper in 
the ACDM model, due to the fact that fewer galaxies tend 
to form on average at low redshifts in such a model. More 
specifically, this increase in the steepness of the SFR results 
from: 

• a similar high redshift SFR level in a ACDM and in a 
SCDM cosmology because the earlier formation of galaxies 
in a ACDM model is compensated by a higher normalization 
of the power spectrum in a SCDM model 

• a lower contribution to the SFR history at late times in 
the ACDM model due to the formation of a smaller number 
of galaxies than in the SCDM model 

This explains why the discrepancy mainly shows up in the 
spiral population and that, in contrast, the contribution of 
elliptical galaxies which form much earlier is very similar in 
both cosmologies (Silk & Devriendt 2000). We argue that the 
amount of star formation occurring in dust-shrouded objects 
cannot be much greater than our model predicts in order 
to avoid overestimating the submm diffuse background (see 
figure 6). Therefore we emphasize that cosmic star formation 
has to decrease at redshifts z > 3.5. Note that this result is 
strengthened if dust-shrouded AGNs are a major contributor 
to the submm emission, even though this does not appear 
to be the case if X-ray characteristics are a reliable AGN 
monitor (e.g. Barger et al. 2001). 



4 CONCLUSIONS 

We have implemented the model of BSS98 within the sim- 
ple semi-analytic model of DG00 in order to describe galaxy 
collisions on a physical basis. We find that such a combina- 
tion: 

• naturally reproduces galaxy counts at various wave- 
lengths and the diffuse background radiation. 

• yields a fair match to existing redshift distributions 
ranging from the optical (CFRS) to far-IR (IRAS at 60 mi- 
crons). Furthermore, the model also predicts redshift distri- 
butions of sources at still larger wavelengths (ISOPHOT at 
175 microns; SCUBA at 850 microns) 

• predicts the contribution of morphological types to 
the previous quantities, with spheroids strongly dominat- 
ing bright galaxy counts only in the sub-millimetre window 
(from 450 microns and longwards). 

• predicts that the cosmic star formation rate history 
is dominated by spiral-like objects until redshift 2 where 
spheroids become the major contributors. 

Our results qualitatively agree with the model of the infrared 
universe presented in Tan, Silk & Balland (1999). We also 
note good agreement of the main conclusions of the present 
work with the collisional starburst model of Somerville et 



al. (2001). This agreement is examplified in the comparison 
between our SFR plot (figure 8 of the present paper) and 
the SFR predicted by their collisional starburst model (their 
figure 9). It is satisfying that these two approaches give con- 
sistent results as it supports the view that galaxy-galaxy 
collisions in the high redshift universe might have played a 
dominant role in triggering star formation. 

Note however that Somerville et al. fail to account 
for the submm counts (Somerville private communication). 
They find too few highly obscured luminous galaxies. Our 
model gives satisfactory agreement both in the restframe 
blue and FIR. This is because we have the freedom, and in- 
deed the motivation from the model, to cumulate the effects 
of collisions. In contrast, semi-analytic hierarchical models 
form spheroids over a Hubble time at the redshift at which 
galaxy mergers occur. This may not be enough if the dust 
is dispersed by the ensuing starburst. Introduction of an 
ad hoc delay between multiple starburst trigerred by differ- 
ent orbital stages during the pre-merger phase might solve 
this. Alternatively one might appeal to additional sources 
of IR luminosity such as the one that might be associated 
with AGN formation. All indications however suggest that 
AGNs do not play a major role in accounting for most of the 
observed submm sources, as evidenced from the CHANDRA 
deep fields. In our case we have a physical model that works 
in the right direction and yields a self-consistent panchro- 
matic explanation of the observed universe both locally and 
at high redshift. We are aware that what we called spheroids 
at high redshifts but are in the process of a violent burst of 
star formation will not have the morphologies of a classical 
early type galaxy, as it will take some time for these objects 
to relax. We defer a more detailed study of the evolution of 
relaxed galaxies to a companion paper, but note that it is 
likely that this fraction will be fairly small at z > 3 since 
there will be very little time elapsed since the collision. 

However, in spite of this issue, an important conse- 
quence of our model, is that the FIR counts and luminos- 
ity density are expected to continue to rise with redshift to 
z ~ 3. Hence we predict that a future mission with adequate 
resolution to avoid source confusion, which sadly may not 
be the case for ASTROF or SIRTF, will provide an important 
test of the galaxy collision model. 
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